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Abstract: The current research around the Fibonacci's and Lucas' sequence evidences the 
scientific vigor of both mathematical models that continue to inspire and provide numerous 
specializations and generalizations, especially from the sixthies. One of the current of research and 
investigations around the Generalized Sequence of Lucas, involves it's polinomial representations. 
Therefore, with the introduction of one or two variables, we begin to discuss the family of the 
Bivariate Lucas Polynomias (BLP) and the Bivariate Fibonacci Polynomials (BFP). On the other 
hand, since it's representation requires enormous employment of a large algebraic notational 
system, we explore some particular properties in order to convince the reader about an inductive 
reasoning that produces a meaning and produces an environment of scientific and historical 
investigation supported by the technology. Finally, throughout the work we bring several figures 
that represent some examples of commands and algebraic operations with the CAS Maple that 
allow to compare properties of the Luc as'polynomials, taking as a reference the classic of 
Fibonacci's model that still serves as inspiration for several current studies in Mathematics. 
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1. Introduction 


Undoubtedly, the Fibonacci sequence preserves a character of interest and, at the same time, of 
mystery, around the numerical properties of a sequence that is originated from a problem related to the 
infinite reproduction of pairs of rabbits. On the other hand, in several books of History of Mathematics 
in Brazil and in the other countries (Eves, 1969; Gullberg, 1997; Herz, 1998; Huntley, 1970; Vajda, 
1989), we appreciate a naive that usually emphasizes eminently basic and trivial properties related to 
this sequence. This type of approach can provide a narrow and incongruent understanding of the 
Fibonacci sequence, mainly about it's current evolutionary stage. 


On the other hand, from the work of the mathematicians Francjois Edouard Anatole Lucas (1842 - 
1891) and Gabriel Lame (1795 - 1870), we observe a progressive return by the interest of the study of 
numerical sequences and their properties. Thus, we highlight the following set of the numerical 


sequences: {l, 1,2,3,5,8,13,21,34,55,89,144,233,377,610,987,1597,2584,4181,6765,10946, 

{1,3,4,7,11,1 8 , 29,47,76,123,199,322,521, m3, 1364,2.207,3571,5778,9349,15127,. .., L,,...} 
;{1,2,5,12,29,70,168,408,,..„P, L,= L, +L„_,; P, tl =2P n _ t +L„_ 2 . n>\. 


We observed that the initial values are indicated by f x = 1, f 2 = 1, L, =, L 2 = 3, Pj = 1, P, = 2 . The 

Fibonacci polynomials were first studied in 1883 by Belgian mathematician Eugene Charles Catalan 
(1814 - 1894) and German mathematician Ernest Erich Jacobsthal (1881 - 1965). Thus, Catalan 
defined the following family of Fibonacci polynomial functions as follows. 

Definition 1: We will call the Fibonacci Polynomial Sequence - SPF, the set of polynomial functions 
described by the recurrence relation f ] (A) = 1, f 2 (A) = x, f n (x) = x ■ f n _ x (x) + f n2 (A), n > 1. 

Undoubtedly, through the previous definition, we can perceive the generalization process that initially 
occurred with the Fibonacci sequence and, gradually, after a few decades, began to be registered in 
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several other generalized sequence models as well. Thus, in the following section, we will study the 
recursive relations arising from the introduction of two real variables, even when we consider the 
introduction of an imaginary unit. Such mathematical definitions are relatively recent and reflect an 
evolutionary mathematical and epistemological process of the Fibonacci and Lucas" sequence. 

2. The Bivariate (Complex) Fibonacci and (Complex) Lucas Polynomials 

The scientific interest in returning study and research around the Fibonacci and Lucas sequences can 
be derived from the inauguration and publication of the periodical entitled The Fibonacci Quarterly, in 
1963 (Gullberg, 1997, p. 243; Posamentier & Lehmann, 2007, p. 28; Stakhov, 2009, p. 130). From this 
fact, there was a greater publicity and dissemination of information and scientific data related to a 
diversity of repercussions and progressive specializations in the research subfields of the Lucas and 
Fibonacci model (Brousseau, 1965, 1967; Gould, 1981; Honsberger, 1985; Livio, 2002; Vajda, 1989; 
Walser, 2001). 

As we mentioned in the previous section, one of the forms of specialization of the Fibonacci model 
occurred from the introduction of a kind of parameter or variable and that gave rise to the definition of 
the Fibonacci polynomial sequence (see definion 1). Similarly, from the seventies, a similar kind of 
generalization of the Lucas" model became better know. 

For example, Asci & Gurel (2012, p. 1) mention that the greater interest in the study of the Lucas 
polynomials in 1970 occurred with the work of Bicknell (1970) and Hoggat & Bicknell (1973a; 
1973b). In addition, for it"s formulation, we recorded the following recurrence. 

Definition 2: We will call the Lucas Polynomial Sequence - SPL, the set of polynomial functions 
described by the recurrence relation 

L^(x) = 2, L, (x) = x,L ti (x) = x ■ L n i (x) + L n 2 (x) . (Bicknell, 1970). 

From a tradition of works (Hoggat & Long, 1974; Swamy, 1968; Web & Parberry, 1969) that 
discussed properties of the bivariate polynomials, we still have the following definition. 

Definition 3: The Generalized Bivariate Fibonacci Polynomial (GBFP) may be defined as 
H 0 (x, y) = a 0 , H, (x, y)=a 1 , H„ +1 (x, y) = x ■ H n _ x (x, y) + y ■ H n 2 (x, y), n > 1 . We assume 

A, y ^ 0, x 1 + 4y ^ 0. (Catalini, 2004a; 2004b). 


From a particular cases, taking a 0 = 0, a x =1 and a 0 = 2,a l = x we can define. 


Definition 

4: 

The 

Bivariate 

Fibonacci 

Polynomial 

(BFP) may be 

defined 

F 0 (x, y) = 

0,Fj(x, y) 

= 1, F +1 (x, y) = x • F n _ 

i (x, y) + y • F, 

_ 2 (x,y),n>l. 


Definition 

5: 

The 

Bivariate 

Lucas 

Polynomial 

(BLP) may be 

defined 


A)(■*> y) = 2 > L 1 (*> y) = L„ + 1 (x, y) = x ■ L n _ x (x, y) + y • F„_ 2 (x, y ), n > 1 . 

We introduce the unity imaginary ‘i’ in order to define the following sequences. 

Definition 6: We will call the Bivariate Complex Fibonacci Polynomial Sequence (BCFP) described 
by the recurrrente relation F 0 (x, y) = 0, F, (x, y) = 1, F n+l (x, y) = ix ■ F n (x, y) + y • F n _ x (x, y), n> 1. 

In the similar way, Asci & Gurel (2012) bring a last definition 

Definition 7: We will call the Bivariate Complex Lucas Polynomial Sequence (BCLP) described by 
the recurrrente relation L fl (x, y) = 2, L[ (x, y) = ix, L ii+l (x, y) = ix ■ L n (x, y) + y • L n , (x, y), n> 1. 
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From the previous mathematical definitions, we can easily determine the polynomial forms (BFP, 
BLP, BCFP, BCLP) we can see in the figure below. Clearly, we note the presence of the imaginary 
unity of algebraic expressions that tend to become increasingly complex, as the order of the 
corresponding element increases. Thus, similarly to the historical elements that acted as 
epistemological obstacles to the establishment of an explicit formula of the Fibonacci and Lucas terms 
(Spickerman, 1982; Waddill & Sacks, 1967; Witford, 1977), we may express a similar interest for the 
determination of an explicit formula of the terms presented in the bivariate polynomial sequence. 


On the other hand, although we find the effort of some authors to formulate certain algebraic 
representations that may facilitate the determination of particular elements F n (x, _y)„ c/;V and 

L n (x, y)„ c /,v of any order present in the sequence, we also find the high operational cost conditioned 

by the use of high order matrix representations. Consequently, when we explore the computational 
model we can anticipate and test certain properties and the validity of certain mathematical theorems. 
In this case, Asci and Gurel (2012) comment the first elements present in the Bivariate Complex 
Fibonacci and Lucas" family, and then in Figure 2, we can observe the proposition of a matrix, of 
order ’n’, that allows determining explicitly any of it's elements. 


n 
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and 
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Figure 1. Ascii & Gurgel (2012) list the first terms of the Polynomial Bivariate Fibonacci and Lucas sequence. 
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Figure 2. Ascii & Gurgel (2012) indicate the matrices whose determinants produce the elements of both 

sequences 


In a precise way, on the left side, the determinant of the matrix allows to ascertain any element of the 
type F n (x, y ) = det(D ;i (x, y)) , for n> 1, while on the right side, through the determinant of the 

indicated matrix, we can make explicit from the formula L n (x, y) = det(H n _j (x, y)), n> 1 (see 


Volume 9 Number 4, 2016 
































74 


Francisco Regis Vieira Alves, Paula Maria Machado Cruz Catarino 


Figure 2). Certainly, when we deal with matrices of fourth, fifth, and sixth and N-th orders, we will 
have serious operational difficulties and, in view of this, the necessity of using the computational 
model is unquestionable. Before, however, we will see in the next section some elemental properties 
of divisibility related to the Lucas's sequence, that we have indicated by j L n . 


4 + i 4 
v 4 4-i j 


3. Some preliminarly properties 

Jeffery & Pereira (2014) take the particular matrix A = 
induction, it will be easy to determine that A Q" — 
hand, we still observe A = 

following identities A 2 =51, A- Q = Q ■ A = 

Lemma 1: Let k,n e IN , with k odd, then L \L, . 

’ ’ ’ n k-n 


1 2 
v2 -4 


and Q = 


4 4 
V 1 o y 


Thus, by 


= Q" ■ A, for n > \. On the other 
= 2Q — I. Moreover, we have the 

4 n (h 4 A 


4 2 ^ 

= 2 

4 n 
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(4 

+2 -4 


4 <4 
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(commutative property). 


Proof. Preliminary, we take AQ' = 


expression, we obtain deX(AQ") = det 


f J J 'N 

+n+l n 

v 4 4-i j 
4+i 4 
v 4 4-i j 


and, taking the determinant of the same 


= LT ,—L 2 . But we will also have to 

n +1 n —1 n 


det(A0”) = det(A)det(2") — — 5-det(Q)---det(0) =5-(—l)" +l . Now, if we consider that, for a 

particular value 'n Q ' we have the propertie 5 \ and, but in view of 5 \ L t] L^ , — L n 2 , we could 

conclude that 5\ L +| L„ n i 5 \ L nn { . Finally, we must verify that 5\L j =L +L l . However, 

this can not occur. Thus, no element present in the Lucas sequence is divisible by 5. On the other 
hand, since k is odd, we can write k = 2m +1. Now, we observe that the matrix 


AQ" = 


^ l n+\ ^n 

v 4 4-i j 


is a diagonal matrix mod L n . Morover, we get the matrix 


AQ' ■ • • AQ' = ( AQ" ) k is also a diagonal matrix mod L n . In addition, we still observe that holds 
(AQ") k =(AQ") 2m+l =A 2m+1 Q" k =(A 2 )'"A-Q" k =(5 m I)A-Q nk =5 m AQ" k . From this, we 


obtain the identity 5" 


~ J nk +1 ^nk 


v 4* 4t-i j 


= 5"' ■ AQ" k =(AQ") k is also a diagonal matrix modL . 
Since L is not divisible by 5, we must have L \L , , with k odd. 

n J 7 n nk ’ 

Lemma 2: Let k , n e IN , with k even, then L \L, . 

’ 5 ’ n k-n 

Proof. In a similar way, we record (AQ" f = (AQ" ) 2m = A 2m Q nk = (A 2 )"' ■ Q" k = (5'"I ) • Q" k = 5'" ■ Q" k . 


From this, we observe 5 m ■ 


f j L ^ 

^nk+l *-•nk 


V nk 


4/t-l 


= 5"' • Q" k = (AQ n ) k . Again, we repeat the argument that, 


we observe that the matrix AQ" = 


L +. L 

72+1 77 

v 4 4-i j 


is a diagonal matrix mod L n . Morover, we get the 
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matrix 5'" ■ Q" k = (AQ") k is also a diagonal matrix modL ;( . Follow the result, Since L n is not 
divisible by 5, we must have L n \ L nk , with k even. 


m n 


Theorem 1: Let m ,n & IN and let d = gcd(/«, n). If the numbers —,— are both odd, then 

d d 

gcd (L m ,L n ) = L d . 


Proof. From the lemma 1 and 2, by definition we know d \m,d\n and, consequently, we have 
L d \ L n and L d \ L m . Now, we show that gcd(L m ,L n )\ h d . If we consider the particular case 
gcd(L m , L n ) = 1 \ L^, it's done. So we suppose gcd(L m , L n ) = d > 1. We know some basic properties 
that we can get the intergers numbers a,b such that am + bn = gcd(m, n) (Koshy, 2007; Tattersall, 
2005). We still observe that if ‘a’ is odd, the other must be even (b is even), or reciprocally, in view of 
the properties of the G.C.D.. So, we take the matrices AQ" l ,AQ" are both diagonal matrices 


mod[gcd( L m , L n )]. On the other hand, we consider the matrices (AQ" 1 j , (AQ ") and the 
matricial product indicated ( AQ'") -( AQ" j = A“Q‘" n ■ A b Q b " = A“ . ^Q am+bn _ We record that 
2p + l,b = 2q and, follows that ( AQ’"J ■(AQ")’ =(A 2 ) p ■ A-(A 2 ) q Q am+bn =(5I) P A(5I) q )-Q d 


a = 


= 5 


p+q 


■ {A ■ Q d ). Finally, we still observe that ( AQ m J ■ [aQ"J = 5 p+q ■ (A • Q d ) = 5 p+q 


Lj+i L d 
y L d L d-\j 


But, the matrices AQ m ,AQ n are both diagonal matrices mod[gcd(L m ,L n )] and, the same propertie 
we have for the both matrices ( AQ,( AQ" j . Hence, |gcdl L fu ,L n )\ divides and 


gcd (L m , L „) = d. 


In this section we study properties involving the divisibility of the numbers present in the Lucas 
sequence. Naturally, similar properties are widely known related to the Fibonacci sequence. In 
addition, from a specialized literature, we also know that the properties of divisibility are preserved, 
since we deal with polynomial functions in one or in two variables. Surprisingly, we shall see that the 
same regularity is not expected for the set Fibonacci numbers and the Lucas polynomials, (see fig. 3). 
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FI02 = 2 * 1597 * 3571 * 6376021 ** 919 * 3469 

FI04 = 3 * 7 =*= 233 * 521 * 90481 ** 103 * 102193207 
FI 05 = 2*5*13*61* 421 * 141961 ** 8288823481 
F106 = 953 * 55945741 ** 119218851371 

FI08 = 2 4 * 3 4 * 17 * 19 * 53 * 107 * 109 * 5779 ** 11128427 

FI 10 = 5 * ll 2 * 89 * 199 * 661 * 474541 ** 331 * 39161 

FI 12 = 3 * 7 2 * 13 * 29 * 47 * 281 * 14503 ** 10745088481 
FI 14 = 2 3 * 37 * 113 * 797 * 9349 * 54833 ** 229 * 95419 
FI 16 = 3 * 59 * 19489 * 514229 ** 347 * 1270083883 

FI 17 = 2 * 17 * 233 * 135721 ** 29717 * 39589685693 

FI 18 = 353 * 2710260697 ** 709 * 8969 * 336419 

FI20 = 2 5 * 3 2 * 5 * 7 * 11 * 23 * 31 * 41 * 61 * 2161 * 2521 ** 241 * 20641 

FI26 = 2 3 * 13 * 17 * 19 * 29 * 211 * 421 * 35239681 ** 1009 * 31249 

FI 28 = 3 * 7 * 47 * 1087 * 2207 * 4481 ** 127 * 186812208641 

FI 29 = 2 * 433494437 ** 257 * 5417 * 8513 * 39639893 

FI 30 = 5 * 11 * 233 * 521 * 14736206161 ** 131 * 2081 * 24571 

FI 32 = 2 4 * 3 2 * 43 * 89 * 199 * 307 * 9901 * 19801 ** 261399601 

FI 34 = 269 * 116849 * 1429913 ** 4021 * 24994118449 

FI 38 = 2 3 * 137 * 139 * 461 * 829 * 18077 * 28657 ** 691 * 1485571 

FI 40 = 3 * 5 * 11 * 13 * 29 * 41 * 71 * 281 * 911 * r41961 ** 12317523121 


Figure 3. Dresel & Daykin (1965) studied the factorization of the Fibonacci's numbers for n>100. 


However, before we begin the next section, we recall the following property numerical behavior 
f l9 = 4181 = 37 • 113. Vemner & Hoggatt (1974, p. 113) mention that, in the Fibonacci sequence, the 

condition that ‘n’ be a prime is necessary but not sufficient for the primarity of the corresponding 
value. Moreover, the research about the primarity test relatively the both sequences has been 
developed for decades (Brillart; Montgomery & Silverman, 1988; Daykin & Dresel, 1970; Dubner & 
Keller, 1999). For example, Daykin & Dresel (1970, p. 30) observe that L,., =151549-11899937029, 

while 73 is a prime. In Figure 3, we can identify other cases wth the support of the technology. On the 
other hand, we will see that some desired properties can be verified in the case of the polynomial 
functions discussed in this work (Daykin & Dresel, 1970). (see Figure 3). 

However, we record the definition 1, and the same authors declare that f (x) is irreductible over the 
ring Z[x\ if only if ‘p’ is a prime. In addition, we have the following property over the ring Z[x, y ]. 

Lemma: If p(x, y),q(x, y),r(x, y) e Z[x, y\, p(x, y) \ q(x, y) ■ r(x, y) and p(x, y) is irreductible, 
then p(x,y)\q(x,y) or p(x, y) \ r(x, y). 

This will be discuss in the next section with the Maple's help. 

4. Some properties of the Bivariate Lucas Polynomial Sequence with the Maple 's help 

We can mention pioneering works that discussed some divisibility criteria involving polynomial 
functions obtained by recurrence, present in the Fibonacci and Lucas sequence, in one or several 
variables (Bicknell, 1970; Hoggatt & Bicknell, 1973a; 1973b; Hoggatt & Long, 1974; Webb & 
Parberry, 1969). In particular, we note the results discussed Hoggatt & Long (1974). In fact, we 
enunciate the following results, in accordance the formal definition 
F 0 (x, y) = 0, F, (x, y) = 1, F n+] (x, y) = x ■ F n (x, y ) + y ■ F n _ x (x, y), n > 1. With origin in these works 
of the seventies, we will announce some important results. 
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Theorem 2: For m>0,n>0 we have F m+n+l (x,y) = F m+I (x, y)F n+l (x, y) + y- FJx,y)F n (x, y ). 
(Hoggatt & Long, 1974, p. 114). 

Proof. We observe that FJx,y) = F 0+l Jx,y) = x-F i (x,y) + yFJx,y) = F 0 Jx,y)-F M {x,y) + y-FJx,y)F l (x,y). 
Moreover, for n — 2 we write F 3 (x, y) = F 1+1+1 (x, y) = x ■ FJx, y) + y ■ F x (x, y) = x ■ x + y • 1 • 1 = 
= F 2 (x,y)F 2 (x,y) + y-FJx,y)FJx,y) = F u Jx,y)F l+ Jx,y) + y-FJx,y)FJx,y). For a fixed 
integer m> 0 we proceed for induction for ‘n’ and we assume that the propertie is true for 
F m+n+ Jx,y) = F m+ Jx,y)F n+ Jx,y) + y- FJx,y)F n (x,y). In the next step, we will consider the 

element F m+n+ Jx,y) = x-F m+n Jx,y ) + y-F m+ Jx,y) = x-(F n Jx,y)F n Jx,y) + y- FJx,y)FJx,y)) + y-F m Jx,y). 
But, we observe that F m+ Jx, y) = F m+(n _ ij+i (x, y) = F m+X (x, y)FJx, y) + y- FJx, y)F n _Jx, y). 
Now, we take xF m+x (x, y)F +l (x, y) + xy- FJx, y)FJx, y) + y • F n+I (x, y)FJx, y) + y 2 - FJx, y)F n Jx, y) = 
= l*F m+ 1 (x, y)F n+1 (x, y) + y- F m+1 (x, y)F n (x, y)] + [x y- F m (x, y)F n (x, y) + y 2 • F m (x, y)F n _ x (x, y)] = 
= F m J x , y)[x-F n+l (x, y) + yFJx, y)] + y- FJx, y)[xFJx, y) + y- F n Jx, y)] = F m Jx, y)F (n+V) Jx, y) + 
+ y- FJx,y)F n Jx,y). Finally, we obtained that: F m+n+1 (x, y) = F m+] (x, y)F in+V)+] (x, y) + y- FJx, y)F n+] (x, y), 
for every m > 0, n > 0. 

Lemma: For n> 0 , then gcd(y, FJx, y)) = 1. (Hoggatt & Long, 1974, p. 114) 

Proof. The assertion is clearly true for n = 1gcd(y, FJx, y)) = gcd(y,l) = 1. Assume that it is true 
for any fixed k> 1. Then, since F k+ Jx,y) = x- F k (x,y) + y ■ F k _Jx,y) . If the condition does not 
hold, we could take gcd(y, F k+ Jx, y)) = d(x, y). But, en virtue the definion, we have 
d(x,y)\F k+ Jx,y), d(x,y)\y. Therefore, we will have d(x,y )\F k+ Jx,y)-y-F k Jx,y) = x-FJx,y). 
However, it can not be d(x,y )\FJx,y), since gcd(y, FJx, y)) = 1. Therefore, we obtain that 
d(x,y)\x but, this can not occur unless that gcd(y, F k+ Jx, y)) = 1. 

Theorem 3: For n > 0 , then gcd(F (x, y), F n+ Jx, y)) = 1. (Hoggatt & Long, 1974, p. 116) 

Proof. Again, the result is trivially true for n = 0, n = 1 since that gcd( FJx, y), FJx, y)) = gcd(0,1) 
and gcd( F { (x, y), F ? (x, y)) = gcd(l, x) = 1. We assume that is true for n = k — 1, where ‘k’ is a fixed 
interger k> 2 and we assume that gcd (F k (x,y),F k+ Jx,y)) = d(x,y). Since we know 
F k+ 1 (x, y) = x • F k (x, y) + y • F k l (x, y) . Again, by definition of the g.c.d., we will have the propertie 
d (x, y) \ F k+1 (x, y) - x • F k (x, y) = y • F k , (x, y) . Thus, we will have d (x, y) \ y ■ F k , (x, y). But, we 
know that gcd( F k (x, y), F k , (x, y)) = I and, in this way, F k _Jx,y) is not divisible by d(x, y). 
Consequently, d(x, y) \ y that is a irreductible polynomial in the variable ‘y’. So, it's a contradiction. 

Theorem4: For m> 2, then F m (x, y) \ F (x, y) <=> m\ n . (Hoggatt & Long, 1974, p. 116) 

Proof. (=>) Again, we proceed by induction. Preliminarily, we observe that FJx, y) \ F k Jx, y) is 
true for k = \:.F m (x, y) \ F Un (x, y) . We proceed by induction, for a fixed integer k > 1, that is, we 
know that FJx,y)\F k Jx,y), for k> 1. In the next step, we will see if the following division 
occurs FJx,y)\F {kt] )m (x, y). On the other hand, en virtue the identity of the 
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theorem F m+n+l (x, y) = F m+1 (x, y)F n+l (x, y) + y• F m (x, y)F n (x, y), we write 

F hn (f y)F m+l (x, v) + y- F^, (x, >’)F (x, y). But, since we have (by hypothesis) F m (x, y) \ F k , m (x, y) and 
F m (x,y)\F m (x,y) immediately we obtain FJx,y)\ F km (x,y)■ F m+I (x, y) + y-F fen _,(x,y)• F m (x, y). 
Finally, we find that F n (x, y)\F a+|)m (x, y) . This propertie clearly implies that F m (x, y)\F k m (x, y) , 
for > 1. Thus, we still have m\n=> F m (x,y)\F ii (x,y). (<=)Now, for m > 2 we suppose 

F m (x, y)\F n (x, y) and we must obtain that m\n. On the other hand, if we assume that ‘n’ is not 
divisible by ‘nT, by means of the division algorithm, exist intergers ‘q’ and ‘r’ with the condition 
n = m-q + r,0 < r < m. Again, by the previous theorem, we take the 

formula F n (x, y) = F m . q+r (x, y) = F m . q+(r _ l)+1 (x, y) = F m . q+l (x, y)F r (x, y) + y- F mq (x, y)F r _ l (x, y). 
Finally, we observe that F m (x, y)\ F mq (x, y) by the first part of the proof and F m (x,y)\ F n (x,y) . 
Consequently, we still have FJx, y) \ F (x, y) -y- F nvq (x, y) ■ F ,(x, y) = F ntq+[ (x, y) • F (x, y). 
But, since we know that gcd(F m? (x, y),F +1 (x, y))=l and, the only possibility is 
F m (x,y)\ F r (x,y) but it's impossible, since we have the condition 0 < r < m , that is, the term 
F (x, y) is a lower degree than F m (x, y) . Therefore, r = 0 and m \ n and the proof is complete. 

Theorem 5: For m>0,n>0 we have gcd(F m (x, y), F n (x, y)) = F gcd(m n) (x, y) . (Hoggatt & Long, 
1974, p. 116) 

Proof. Through elementary properties, we know that exist integers r and s, say, r > 0 and ,v < 0 such 
that gcd(m, n) = r ■ m + s n r ■ m = (m,n)-s-n. Thus, by theorem, we write F rm (x, y) = F (mn) _ sn (x, y) 

= F (mn)+( _ s>n (x, y) = F (mn) (x, y)F ( _ s>n+1 (x, y) + y • F^^Cx, y)F ( _ s>n (x, y) = F (m >fl) (x, y)F_^ +1 (x, y) + 
y * F (m n j_\ (x, y) • F_ v ^(x, y). But, let d(x,y) = gcd(F;,(x, y),F n (x, y)) and, consequently we have 
d(x, y) \ F ; . m (x, y) and d(x, y) \ F w! (x, y). From these properties, we still obtain that 
d(x, y) \ F r . B (x, y) - y ■ ^ H (x, y) • F f .(x, y) = F M (x, y) • F wl (x, y) . Now, if we conclude that 
d(x, y) \ F (m n) (x, y) the proof is complete. However, if occur d(x, y) \ F_ yn+l (x, y), we must observe that 
1 = gcd(d(x, y), F m+1 (x, y)). Otherwise, we would have d Yx, y) = gcd(d(x, y), F vn ,, (x, y )) and, 
we still know 1 = gcd(F_ m (x, y), F_ m+1 (x, y)) and d(x, y) \ F_ s . n (x, y). Consequently, we can get 
that d'(x,y)\F. J(1 (x, y) and d'(x, y) \ F s JJ+1 (x, y) and we must obligatorily have to 
1 = gcd(d(x, y), F_ sn+1 (x, y)). Finally, the property d(x, y) \ F im n) (x, y) is verified. 

We will give some examples below in order to verify the expected behavior of some particular cases, 
according to the theorems we have just demonstrated in detail. 

F 19 (x, y) = x 18 + 17x 16 y + 120x 14 y 2 + 455x 12 y 3 + 1001x 10 y 4 + 1287x s y 5 + 924x 6 y 6 + 330x 4 y 7 
+45x 2 y 8 + y 9 (is irreductible over Z[x,y],p=19) 

F 23 (x, y) = x 22 + 21x 20 y + 190x 18 y 2 + 969x 16 y 3 + 3060x 14 y 4 + 6188x 12 y 5 + 8008x 10 y 6 + 6435x 8 y 7 
+3003x 6 y 8 +715x 4 y 9 +66x 2 y 10 + y n (is irreductible over Z[x,y],p=23) 

F 31 (x, y) = x 30 + 29x 28 y + 378x 2 6y 2 +2925x 24 y 3 +14950x 22 y 4 +53130x 20 y 5 +134596x 18 y 6 + 
245157x 16 y 7 +319770x 14 y 8 +293930x 12 y 9 +184756x 10 y 10 +75582x 8 y" +18564x 6 y 12 + 
2380x 4 y 13 +120x 2 y 14 + y 15 (is irreductible over Z[x,y],p=31) 
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Before finalizing the list of some important results, we will enunciate a last theorem that provides an 
important characterization only for the elements present in the Polynomial Fibonacci's sequence. 

Theorem 6: Let r-r(x,y ) be any polynomial in the variables ‘x’ and ‘y\ If there exists a least 
positive interger ‘m’, such that r(x,y)\F m (x, y) , then r(x, y) \ F (x, y) <=> m\ n . (Hoggatt & Long, 
1974, p. 117). 

Proof. By the theorem 4, we know if m\ n => F m (x, y) \ F (x, y). So, we admit that exists a least 
positive interger ‘m’, such that r(x, y) \ F m (x, y) and, by transitivity, we conclude r(x, y) \ FJx, y). 
Now, we suppose that r(x,y)\F n (x, y) and yet ‘n’ is not divisible by ‘m’. Then, by the Euclidean 
Algorithim, exist integers q,s,0 < s < m and n — m-q + s ■ Again, by the theorem 2, we can write 
F„(x, y) = F nvq+s (x, y) = F m . q+(s _ 1)+l (x, y) = F m . q+l (x, y )F s (x, y) + y ■ F m . q {x, y )F s _ l (x, y) - But, since 
r(x, y) \ F m (x, y), r(x, y) \ F (x, y) and r{x, y) \ F (x, y). Consequently, by the last identity, we still 
have r(x, y) \ F (x, y ) - y - F m , q (x, y)F_, (x, y) = F m , q+1 (x, y) F (x, y) . From this, follows that 
r(x,y)\F mq+1 (x,y)F s (x,y) and we know gcd(F mq (x, y), F m , q+1 (x, y)) = 1. Thus, the only alternative 
is r(x, y) \ F s (x,y), however is a crontradicion, since exists a least positive interger ‘m’, such that 
r(x, y) \ F„(x, y) and 0 < s < m . So, ‘n’ is divisible by ‘nT and the proof is complete. 


Now, we will study some properties of the divisibility and factorization of the polynomial elements 
present in the Lucas sequence. We will see that it does not enjoy a similar algebraic behavior. Before, 
however, let us look at some properties of the matrices. 


From the definition L ;; (x, y) = x ■ L n _ l (x, y) + y ■ L n _ 2 (x, y), n > I we can determine some 
particular initial forms: L 0 (.x, y) = 2,F 1 (x, y) = x,L 2 (x, y) = x 2 +2y , L 3 (x,y) = X 3 + 3xy, 
L 4 (x,y) =x 4 + 4x 2 y+ 2y 2 , L 5 (A, y) = +5.r 3 y + 5xy 2 , ,..etc. Catalani (2002) define the matrix 


A(x, y) = 


' x 
vT 


= A and B(x, y) = 




+ 2 y y 

Ay 2y j 


= B ■ Immediately, we will have AB — 


fx n 

( 2 , 0 7 

x + 2 y y 


x 3 +3xy x 2 +2y^ 

, AB 2 = 

^x 4 +4x 2 y + 2y 2 x 3 +3xy ^ 

v>’ °y 

V TV 2 yj 


v (x 2 + 2y)-y x-y j 


v (x 3 +3xy) • y (x 2 +2y)-y y 


x 5 +5x 3 y+5.xy 2 
v (x 4 +2x 2 y + 2x 2 y + 2y 2 )y 


x 4 +4x 2 y + 2y 2 ^ 
(x 3 + 3xy)y y 


x 6 +6x 4 y+9x 2 y 2 + 2y 3 
Jx 5 + 2x 3 y + 3x 3 y+4.xy 2 + xy 2 )y 


x 5 +5x 3 y+5xy 2 
(x 4 +2x 2 y + 2x 2 y + 2y 2 )y y 


In the figure below, we can analise and conjecture a closed form for the product AB",n > 1. So, from 
the appreciation of some particular products, we can acquire a better understanding about the follow 


7 


identity AB" = 




" +2 " +1 , n > 1. And, we can still work with the inverse matrix. 

14+ry L n-yJ 


r 


On the other hand, we can still get the inverse y j — 


7 


0 

1 


A 


V 


x 

yj 


= A -1 and B 4 (x, y) = 
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2 

X 

x 2 +4_y 

x 2 + 4v 

xy 

x 2 +2y 

x 2 +4y 

x 2 +4y 


= B 1 . From this, we can explore the algebraic behavior of the following 


produts: (AB) , (A5) , (A5)" 3 , (AB) , (A,B)“\ (AB) ,..(A#) - ",n>l. Moreover, with 


the CAS, we can find ( AB ) 1 = B 'A 1 = 


( 

X 

x 2 +2y ^ 


{ L,(x,y) 

L 2 (x,y) \ 

x 2 y+4y 2 

(x 2 y+4y 2 )y 


(x 2 +4y)y 

(x 2 +4y)y 2 

x 2 + 2y 

x 3 + 3xy 


L 2 (x, y) 

L,(x,y) 

l v (x 2 y + 4y 2 ) 

(x 2 y + 4y 2 )y j 


v (x 2 +4y)y 

(x 2 +4y)y 2 J 


(AB) 2 = B 2 A 2 = 


f 

x 2 + 2y 


r Lj(x, y) 

L 2 (x,y) ^ 

x 2 y + 4y 2 

(x 2 y + 4y 2 )y 


(x 2 + 4y) y 

(a 2 + 4y) y 2 

x 2 + 2 y 

x 3 + 3xy 


L 2 (x, y) 

L 3 (a, y) 

V (x 2 y + 4y 2 ) 

(x 2 y + 4y 2 )y J 


l (a 2 + 4y) y 

(x 2 + 4y)y 2 J 


In Figure 4, we can observe the behavior of the product of the matrices indicated earlier and, through 
some preliminary cases, by means of an inductive process, formulate it's general term. Of course, we 
can understand that calculus becomes impractical without the use of technology. 

Now let us look at some properties of divisibility and factorization of some of the polynomial 
functions into two variables, present in the Bivariate Lucas Sequence. In Figure 5, we considered a 
strong propertie that permits determine any term of the BPL. In the left side, we can see fourth order 
matrix and, in the right side, we have considered a fifth order matrix. We calculate that 

det(B 4 (x, y)) = Lj (x, y) = x 3 + 3xy = x • (x 2 + 3y) = L, (x, y) • (x 2 + 3y) and det(B 5 (x, y)) = L 4 (x, y) = x 4 + 4x 2 y + 2y 2 . 

We observed that the polynomial in two variables x 4 +4x 2 y+ 2y 2 is irreductible, while we found 
L } (x,y) = L l (x, y) ■ (x 2 + 3 y) . Moreover, we have det (H 9 (x, y)) = (x, y) = x 8 + 8x 6 y + 20 x 4 y 2 + 16x 2 y 3 + 2y 4 
is a another irreductible polynomial over the field /B[x, y]. On the other hand, in Figure 5, we 
observe det(B l3 (x, y)) = L, 2 (x,y) = x 12 + 12x 10 y + 54x 8 y 2 + 112x 6 y 3 +105x 4 y 4 +36x 2 y 5 +2y 6 and, by a 
command of the CAS, we can write det(# 13 (x,y)) = Z 12 (x,y) = (x 4 +4x 2 y+2y 2 )(x 8 +8x 6 y+20x 4 y 2 +16x 2 y 3 + y 4 ) 
= L 4 (x, y)-Z^(x,y), however L, 2 (x, y) is not divisible by L,(x, y) and 12 is not divisible by 8. 
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Figura 4. With the CAS Maple we can verify the algebraic behavior of the product of the matrices 

In addition, with the use of software, we can also determine the factorization and, therefore, the 
decomposition of irreductible factors of the polynomial terms over the ring Z[x,y]. We observe in the 
list below that, unlike the case of Fibonacci, we will have elements of prime index that admit a 
factorization in irreducible factors ( L 5 (x, y) ,L ll (x, y ), L n (x, y)). 

det (H 6 (x, >’)) = L 5 (x, y) = x 5 + 5x 3 y + 5xy 2 = x(x 4 + 5x 2 y + 5 y 2 ) (reductible over Z[x,y],p=5) 

det (H n (x,y)) = L ll (x,y) = x n +llx 9 y+ 44x 7 y 2 + 77x 5 y 3 +55x 3 y 4 + llxy 5 = L,(x,y)(x 10 + llx 8 y + 44x 6 y 2 + 
77x 4 y 3 +55x 2 y 4 +lly 5 ) (reductible over Z[x,y],p=ll) 

det (// 13 (x, y)) = /^(x, y) = (x 4 +4x 2 y+ 2y 2 )(x 8 +8x 6 y + 20x 4 y 2 +16x 2 y 3 + y 4 ) = L 4 (x, y)(x 8 +8x 6 y + 20x 4 y 2 + 
16x 2 y 3 + y 4 ) L 4 (x, y) \ L j2 (x, y),is reductible over Z[x,y] 

det(77 15 (x, y)) = L, 4 (x, y) = x 14 + 14x 12 y + 77x‘°y 2 + 210x 8 y 3 + 294x 6 y 4 + 196x 4 y 5 + 49x 2 y 6 + 2y 7 
= (x 2 + 2y)(x 12 + 12x 10 y + 53x 8 y 2 + 104x 6 y 3 + 86x 4 y 4 + 24x 2 y 5 + y 6 ) = L 2 (x, y)(x 12 + 12x‘°y + 53x 8 y 2 
+104x 6 y 3 + 86x 4 y 4 + 24x 2 y 5 + y 6 ) is reductible over Z[x,y] 
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det(tf i6 (jc,y)) = Z 15 (jc,)0 = x 15 + 15x 13 y + 90x u y 2 + 275x 9 y 3 + 450x 7 y 4 + 378x 5 y 5 + 140x 3 y 6 + 15xy 7 
= x(x 2 + 3y)(x 4 + 5x 2 y+ 5y 2 )(x 8 + 7x 6 y + 14x 4 y 2 + 8x 2 y 3 + y 4 ) = L 1 (x,y)(x 2 + 3y)(x 4 + 5x 2 y+ 5y 2 ) 
(x 8 +7x 6 y+ 14x 4 y 2 + 8x 2 y 3 + y 4 ) is reductible over Z[x,y] 

det (H l7 (x, y)) = L, 6 (x, y) = x 16 + 16x 14 y + 104x 12 y 2 + 352x 10 y 3 + 660x 8 y 4 + 672x 6 y 5 + 336x 4 y 6 
+64x 2 y 7 + 2y s (is irreductible over Z[x,y]) 

det(// 18 (x, y)) = L[ 7 (x, y) = x 17 +17x 15 y+ 119x 13 y 2 +442x n y 3 +935x 9 y 4 + 1122x 7 y 5 + 714x 5 y 6 
+204x 3 y 7 + 17xy 8 = x(x 16 + 17x 14 y +119x 12 y 2 + 442x 10 y 3 + 935x s y 4 +1122x 6 y 5 + 714x 4 y 6 + 204x 2 y 7 
+17y 8 ) (reductible over Z[x,y],p=17) 

det(// l9 (x, y)) = Lj g (x, y) = x 18 +18x 16 y + 135x 14 y 2 +546x 12 y 3 +1287x 10 y 4 + 1782x s y 5 + 
1386x 6 y 6 + 540x 4 y 7 + 81x 2 y 8 + 2y 9 = (x 2 + 2y)(x 4 + 4x 2 y + y 2 )(x 12 + 12x 10 y + 54x 8 y 2 + 

112x 6 y 3 + 105x 4 y 4 + 36x 2 y 5 + y 6 ) = L 2 (x, y)(x 4 + 4x 2 y + y 2 )(x 12 + 12x 10 y + 54x s y 2 +112x 6 y 3 
+105x 4 y 4 +36x 2 y 5 + y 6 ) is reductible over Z[x,y] 

det (H 20 (x, y)) = L, 9 (x, y) = x 19 + 19x 17 y + 152x 15 y 2 + 665x 13 y 3 + 1729x u y 4 + 2717x 9 y 5 + 2508x 7 y 6 
+1254x s y 7 + 285x 3 y 8 + 19xy 9 = x • (x 18 + 19x 16 y + 152x 14 y 2 + 665x 12 y 3 + 1729x 10 y 4 + 2717x 9 y 5 + 
2508x 6 y 6 +1254x 4 y 7 + 285x 2 y 8 +19y 9 ) = Lj(x, y)-(x 18 +19x 16 y+ 152x 14 y 2 +665x 12 y 3 +1729x'°y 4 + 
2717x 9 y 5 + 2508x 6 y 6 + 1254x 4 y 7 + 285x 2 y 8 + 19y 9 )L,(x, y) \ Z 19 (x, y) is reductible over Z[x,y] 

Now, from the list of decomposition into irreducible factors of the elements of the bivariate 
polynomial sequence, we can understand that: L, (x, y) \ L,(x, y) and L s (x, _y) is not divisible by 

L 2 (x, y) = x 2 +2y. Z^(x,y) = x 8 +8x 6 y+ 20x 4 y 2 +16x 2 y 3 + 2y 4 is a irreductible polynomial and 
is not divisible by L,(x, y) = x 2 +2y or L 4 (x, y) = x 4 +4x 2 y+ 2y 2 . In the same manner, the 
bivariate polynomial L l2 (x, y) = x 12 +12x 10 y+54x 8 y 2 +112x 6 y 3 +105x 4 y 4 +36x 2 y 5 +2y 6 is reducitible 
over the ring //?[x,y], since we obtained det(// B (x,y)) = L 12 (x,y) = L 4 (x,y)(x 8 +8x 6 y + 20x 4 y 2 +16x 2 y 3 + y 4 ), 
that is L 4 (x,y)\L l2 (x,y), while L^ 2 {x,y) is not divisible by L,(x, y) or L 6 (x, y). On the other 
hand, we observe L,(x, y)\Z l5 (x, y) and, however, L, 5 (x, y) is not divisible by L 5 (x, y). 
Moreover, we can conclude that L |6 (x, y) is not divisible by the following 
elements L 2 (x, y), L 4 (x, y) and (x, y). And, the only division propertie is L, (x, y) \ L l6 (x, y). 

According to the result indicated by the software, L |7 (x, y) is irreducitible polynomial over the ring 
IR[x, y]. In addition, we further determined that L |9 (x, _y) , that despite having a prime subscript, 
have the element L,(x, y) as a irreducitible factor. Finally, we conclude that L 20 (x, y) is not 
divisible by L, (x, y) , L 5 (x, y) or L |0 (x, y) . In Figure 5, we visualize some command employed en 
virtue to determine it's decomposition over the ring Z[x, y]. 

Now, in addition to an extensive set of algebraic expressions provided by software that indicate the 
decomposition of polynomial functions into two variables, we can draw some conclusions regarding 
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the character of divisibility and factorization of Lucas' polynomial functions. For example, in Figure 5, 
on the left side, we observe that det(// 9 (x,y)) = I 9 (x,y) = x 8 + 8x 6 y + 20x 4 y 2 + 16x 2 y 3 + 2y 4 is 

irreductible, since the software when using the factor command factor [ /, produces the same algebraic 
expression. But, in the same figure, on the right side, we visualize that the element 
det(// 13 (x, y)) = L[ 2 (x, y) = (x 4 +4x 2 y+2y 2 )(x 8 +8x 6 y + 20x 4 y 2 +16x 2 y 3 + y 4 ) = L 4 (x, y) • (x 8 +8x 6 y + 20x 4 y 2 + 
16x 2 y 3 + y 4 ) = L 4 (x, y) ■ L A (x, y) has two components as irreductible polynomials. So, like we have 
mentioned, we know L 4 (x, y)\Z 10 (x, y) , however L,,(x, y) is not divisible by (x, y). 



Figura 5. We obtain the decomposition of the irreducible factors of the elements of the set of polynomials with CAS Maple 

In the next section, we will address an explicit formula for the polynomial terms present in both 
sequences. Yet, we will enunciate some properties involving the Greatest Common Divisor of 

polynomial functions into two variables, now with the introduction of an imaginary unit i 2 = — 1. In 
this way we can compare the class of the BFP (definition 4) with the class of BCFP (definition 6). We 
will find the regularity and invariance of several properties indicated in the theorems discussed here 
and, conversely, the same regularity cannot be observed in the class of BLP (definition 5) and the 
BCLP (definition 7). In fact, we showed that F m (x, y) \ F (x, y) <=> m\n . But, with the software, we 
found several exemples that are counterexamples. 
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Figura 6. We obtain the decomposition of the irreducible factors of the elements of the set of polynomials with CAS Maple 

4. Some properties of the Bivariate Complex Lucas Polynomial Sequence with the 
Maple's help. 

In the previous section, we discuss some properties of bivariate polynomials in two real variables. We 
now turn to the study of a special class of Bivariate Complex Fibonacci and Lucas Polynomials, 
originating from the introduction of an imaginary unit ‘i’ and inherit a tendency of the works 
interested in the process of complexification of the Fibonacci model (Iakin, 1977; King, 1968; Scott, 
1968; Waddill & Sacks, 1967). With this, we can further discuss the process of complexing said 
recursive sequence. Before, however, we recall the definition presented in the predecessor sections en 
virtue to present our first theorem. 

(a'\x, y)-j3'\x, y)) 

Theorem 2: For n>0 we have F n (x,y ) = --- and L n (x, y) = a'\x,y) + fl"(x,y). 

a(x,y)-fi(x,y) 
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Proof. We the characteristic equation designated by t 2 -ix-t — y = 0 . Consequently, we will have the 


following properties and relations between the roots a(x, y ) = 


ix + 




~ x r( ^ 
-, /3(x, y) = 


x 


2 2 
a(x, y) ■ P(x, y) = -y . Finally, from the recurrence relation, we write: F n+l (x, y) = ix- F n (x, y)+y- F n _ , (x, y) = 


a"(x,y) P"(x,y) 

_.{a n {x,y)-p n {x,y))^ {a n ~\x,y)~ p n ~\x,y)) {a n (x,y)~ p'\x,y)) _ a(x,y) /?(*,y) 

uc I *y uc i v 

a(x,y)~P(x,y) " a(x,y)~ P(x,y) a(x,y)~ P(x,y) ' a(x,y)~P(x,y) 

(a"(x,y)-j3 n (x, y)) _ _ (a n (x,y)fi(x,y)-a(x,y)j3 n (x,y)) (a'\x,y)-fi n (x,y)) _ 

i V* ^ , LPv i 

a(x,y)-j3(x,y) ' (a(x, y)fi(x, y))a(x, y) ~P(x, y) a(x,y)~P(x,y) 

| ( a n (x,y)P(x,y)-a(x,y)P''(x,y )) _ ix-a n (x,y)-ix-p l \x,y)-a'\x,y)P(x,y) + a(x,y)P"(x,y) 

{-y)a{x,y)-p{x,y) a{x,y)~ p{x,y) 

_ ix ■ a(x, y) • a" -1 Q, y) - iv • /?0, y) • /3 n ~\x, y)-a"(x, y)/3{x, y) + yj/TQ, y) _ 

a(x, y) - /3(x, y) 

_ a(x, y)” +1 - y • a(x, y)-/3 (x, y )” +1 + y • 1 (x, y) + y • a ”~ 1 (x, y) - y • J3 nl (x, y) _ 

a(x, y)-j3(x,y ) 


q(x, y)” + -/?(*, y) 
a(x, y)-/?(*, y) 


n +1 


. In the similar way, we can determine that L n (x, y) = a"(x, y) + P"(x, y). 


Corollary 1: For the roots of characteristic equation t 2 — ix-t — y = 0 we have the follow identities: 
a „ {x ^_ L n (x,y) + (a(x, y)~P(x, y))F n (x, y) ^ _ L n (x, y)-(a(x,y)~p{x, y))F n (x, y) 


Proof. From the theorem 1, we take the expression ( a(x, y) - P(x, y)) F n (x, y) = a" (x, y) - p n (x, y) 

\ a n (x,y)~ p" (x, yf=(a(x,y)- P(x,y))F n (x,y) 


and LJx, y) = a" ix, y) + P"(x, y) . We solve the system 


a n (x,y) + p"(x,y) = L n (x,y) 


Corollary 2: For n > 0 we have F n (x, y) • L (i (x, y) = F ln ( x , y). 


Proof. Bu the previous theorem, we note F ln (x, y) = 


a-(x, y) - p ln (x, y) _ (a"(x, y) - p n (x, y)) 


/ 


a(x,y)~ P(x,y) a(x,y)~ P(x,y) 


■a" (x, y) + p" (x, y). Follows the result. 

The next results involve some properties related to the divisibility character of the Bivariate Complex 
Fibonacci Polynomials (BCFP). Thus, with respect to the particular case of the numerical sequence, 
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and especially when we introduce two variables, the properties involving the Greatest Common 
Divisor (G.C.D) are invariable preserved over the ring C[x, y]. However, with the computational 
resource, we will see that the same case or pattern does not occur when we deal with the polynomial 
sequence of Lucas. We can compare with the results established in the preceding sections. 

Theorem 2: Let D n (x,y) nxn be a tridiagonal matrices. In this case, by the assumption D 0 (x, y) = 0, 
then we have det D n (x, y) nxn = F n (x, y), n > 0. (Asci & Gurel, 2012). 

Proof. See Figure 2 (on the left side) and more details can be constated in the Asci & Gurel (2013). 
Lemma: Let's suppose that »tr/c(x, y) = 1, n > 0, then gcd(y, F„ (x, y)) = 1. 

Proof. More details can be constated in the Asci & Gurel (2013). 

Teorema: Suposing that gcd(x, y) = 1, n > 0 , then gcd(F (x, y), F M (x, y)) = 1 . 

Proof. More details can be constated in the Asci & Gurel (2013). 

Teorema : For m> 2, we have m, n> 1 F m (x, y) \ F n (x, y) <=> m\ n. 

Proof. More details can be constated in the Asci & Gurel (2013). 

Corollary: Suposing that gcd(x, y) = 1 e m, n > 1 then gcd(F m (x, y), F (x, y)) = F gcd(m n) (x, y) . 
Proof. More details can be constated in the Asci & Gurel (2013). 

However, when we seek to inspect similar properties in the Lucas sequence, from some particular 
cases anticipated by the software, we may notice that several properties are lost and not verified. In 
fact, we list some cases below relatively to the BCLP. 

det(// l3 (x, >’)) = L i2 (x, y) = x 12 -12x 10 y+54x 8 y 2 -112x 6 y 3 +105x 4 y 4 -36x 2 y 5 + 2/ =(x 4 -4x 2 y + 2y 2 ) 
(x 8 -8x 6 y + 20x 4 y 2 - 16x 2 y 3 + y 4 ) = L 4 (x,y) • (x 8 -8x 6 y + 20x 4 y 2 - 16x 2 y 3 + y 4 ) is reductible over Z[x,y] 


det(// 14 (x, y)) = L, 3 (x, y) = x 13 / -13x' 1 yi + 65 x 9 y 2 i - 156x 7 y 3 / + 182x 5 y 4 / - 9 lx 3 y 5 / + 13xy 6 / = 
ix-(x 12 -13x 10 y + 65x 8 y 2 -156x 6 y 3 + 182x 4 y 4 -91x 2 y 5 +13y 6 ) is reductible over Z[x,y] 

det (H l6 (x, y)) = ^(x, y) = -x 15 i + \5x u yi -90x n y 2 i + 275x 9 y 3 i -450x 7 y 4 z +378x 5 y 5 /-140x 3 y 6 / 
+15xy 7 / = -ix-(x 2 -3 y)-(x s -7x 6 y+ 14x 4 y 2 -8x 2 y 3 + y 4 ) -(x 4 -5x 2 y + 5y 2 ) 

is reductible over Z[x,y] 

det {H ll (x,y)) = L lb {x,y) = x 16 -16x 14 y + 104x 12 y 2 -352x 10 y 3 +660x 8 y 4 -672x 6 y 5 +336x 4 y 6 -64x 2 y 7 +2y 8 

is irreductible over Z[x,y] 

det(77 lg (x, y)) = L l7 (x, y) = x 17 / - 17x 15 yz +119x 13 y 2 / - 442x n y 3 z + 935x 9 y 4 / -1122x 7 y 5 / + 714x s y 6 / 
-204x 3 y 7 / + 17xy 8 / = /x-(x 16 -17x 14 y+ 119x 12 y 2 -442x 10 y 3 +935x 8 y 4 -1122x 6 y 5 +714x 4 y 6 - 
204x 2 y 7 +17y 8 ) is reductible over Z[x,y] 
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det(// I9 (jc, y)) = L i8 (x, y ) = -x 18 + 18jc 16 y - 135x 14 y 2 + 546x 12 y 3 -1287jc 10 y 4 + 1782x s y 5 -1386// + 540x 4 y 7 
-81x 2 y s +2y 9 = -(x 2 -2y)-(x 4 -4x 2 y + y 2 )-(x 12 -12x 10 y + 54x 8 y 2 -112x 6 y 3 + 105x 4 y 4 -36x 2 y 5 + y 6 ) = 

= Lj(x, y)-L 4 (x, y)-(x 12 -12x 10 y + 54x 8 y 2 -112x 6 y 3 + 105x 4 y 4 -36x 2 y 5 + y 6 ) 
is reductible over Z[x,y],p=18 

det (H 20 (x, y)) = L l9 (x, y) - -x 19 / + 19x 17 y/ - 152x 15 y 2 / + 665x 13 y 3 / -1729x"y 4 / + 2717x 9 y 5 / 
-2508x 7 y 6 / + 1254x 5 y 7 /-285x 3 y 8 / + 19xy 9 / = -/x-(x 18 -19x 16 y+ 152x 14 y 2 -665x 12 y 3 + 1729x 10 y 4 
-2717x 8 y 5 + 2508x 6 y 6 - 1254x 4 y 7 + 285x 2 y 8 -19y 9 )(is reductible over Z[x,y],p=19) 
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Figura 7. Factoring of terms of the bivariate sequence of Lucas with the CAS Maple (produced by the authors). 


Now, in summary, we show that Lj 3 (x, y) is reductible over □ [x, y]. L l4 (x, y) is divisible by 
L 2 (x, y) = ix, but is not divisible for Lj(x,y). L, 6 (x, y) is divisible by L 2 (x, y) = ix, but is not 
divisible by L 4 (x, y) = -x 3 / + 2xy/. 


On the other hand, L, 6 (x,y) = x 16 -16x 14 y + 104x 12 y 2 -352x 10 y 3 +660x 8 y 4 -672x 6 y 5 +336x 4 y 6 -64x 2 y 7 +2y 8 
is a irreductible polynomial, while Z 18 (x, y) is not divisible by L } (x,y),L^(x,y) and is divisible 
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by L 2 (x, y) = — x +2y and L 4 (x, y) = x 4 — 4x 2 y + y 2 . Finally, L I9 (jc, y) has two irreducible factors 

and none of them belongs to Lucas Polynomial' sequence. Thus, similar to the case of the Lucas' 
sequence we saw in the last section, when we deal with the introduction of the imaginary unit, we 
realize that the same properties that we study with the Fibonacci sequence are lost. 
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Figura 8. Determination of an element of nineteen order in the BLP Lj 9 (a, y) . (produced by the authors). 


4. An historical investigation with the Maple's help 

In some of our works (Alves, 2017; 2016a; 2016b; 2016c; 2015) that seek to emphasize a perspective 
that points out the importance of the History of Mathematics, with a view to an understanding of the 
process of emergence, evolution and systematization of ideas, models, abstract properties, theorems 
and, above all, formal definitions. We have taken a position that Mathematics envolves as a process of 
debugging, synthesizing and testing the adequacy and robustness of certain formal mathematical 
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definitions. Thus, we assume that for an evolution of the evolutionary process in Mathematics, an 
understanding about the process that contributes to such demarcation is essential. 

On the other hand, certain definitions determine a charged notational representation system that 
sometimes prevents an understanding and the establishment of necessary conceptual relations in order 
to identify the invariant elements in each case. With this function, CAS Maple allowed us to explore 
several particular cases, verifying and proving properties originating from an inductive model or 
reasoning. In addition, CAS Maple makes it possible to perform calculations that are impractical when 
we disregard the technology. 

Regarding the Maple's use, we highlight the elements: (i) The software enables verifications of 
particular cases and properties related to the Bivariate Fibonacci Polynomials and the Bivariante Lucas 
Polynomials; (ii) The software allows the verification properties provided by classical theorems 
related to the Bivariate (Complex) Fibonacci Polynomials and Bivariante (Complex) Lucas 
Polynomials, especially the most recently discussed in the literature; (iii) The software allows the 
description of a lot of special particular cases conditioned by newly formulated mathematical 
definitions; (iv) The software enables verification of properties related to a larger set of integer 
subscripts indicated in the scientific articles; (v) the software enables verification of a large number of 
individual cases in order to test mathematical conjectures about the BFP and BLP; (vi) the software 
allows the correction of mathematical formulas in order to provide a precise description. 

Before concluding the current section, we bring a last formal definition that expresses the current 
character of the discussion of the subject addressed to the course of this work. In the last definion we 
constate the Fibonacci' model described by the complex variable z = x + iy . 

Definition 4: Let a Generalized Polynomial {F n (a, z)} in the variables ‘a’ and ‘z’, is designed for 
the recurrence relation F n (a, z) = a- z- F n x (a, z) + a 2 - F n _ 2 {a, z), n > 2 , with initial 
conditions F 0 (a, z) = 0, F y (a, z) = 1. (TASKOPRU & ALTINTAS, 2015). 

In the particular case for the value a = l,/ 2 (z) = Z / 3 (z) = Z“ + l, / 4 (z) = Z +2z = z(z" + 2) , 

/ 5 (z) = z 4 +3z 3 +1, / 6 (z) = z 5 +3z 4 + z 3 +3z = z(z 2 +3)(z 2 +1), f 7 (z) = z +5z 4 + 6z 2 + l, 
/ 8 (z) = z 1 + 6z 5 + 10z 3 + 4z = z(z 2 + 2)(z 4 + 4z 2 + 2), / 9 (z) = z 8 +7z 6 +15z 4 +10z 2 +1 = (z 2 +l)(z 6 +6z 4 +9z 2 +1), 
f 10 (z) = z 9 + 8z 7 + 21z 5 + 20z 3 + 5z = z(z 4 + 3z 2 +1 )(z 4 + 5z 2 +5), f n (z) = z'° + 9z 8 + 28z 6 
+35z 4 + 15r +1, f l2 (z) = z 11 +10z 9 +36z 7 +56z 5 + 35/ +6z = z(z 2 +3)(z 2 +2)(z 2 +l)(z 4 + 4z 2 +1), 
/ 13 (z) = z l2 +1 lz 10 + 45z 8 + 84z 6 + 70z 4 + 21z 2 +1. / 14 (z) = z 13 +12z u +55z 9 +120z 7 +126z 5 +56z 3 +lz = 

z(z 6 + 5 z 4 + 6z 2 + l) (z 6 + 7z 4 + 14Z 2 + 7) , f l5 (z) = z 4 +13z 12 +66z 10 +165z 8 +210z 6 +126z 4 +28z 2 +1 = 

(z 2 + l) (z 4 + 3z 2 + l) (z 8 + 9z 6 + 26z 4 + 24z 2 + l) ; / 16 (z) = z' 5 +14z° +78Z 11 +220z 9 +330z ? +252/ +84/ +8z = 
z(z 2 +2) (z 4 +4z 2 +2) (z 8 +8z 6 +20z 4 + 16z 2 +2) ? / 17 (z) = Z 16 + 15z' 4 + 91z' 2 + 286z 10 

+495z 8 + 462z 6 + 210z 4 + 36z 2 +1, / 18 (z) = z 17 + 16z 15 + 105z 13 + 364z u + 715/ + 792z 7 + 462z 5 + 120z 3 + 9z 
z(z 2 + 3)(z 2 + l)(z 6 + 6z 4 + 9 z 2 + l)(z 6 + 6z 4 + 9 z 2 + 3) ? / 19 (z) = z 18 +17z 16 +120z* 4 + 455z 3i + 
lOOlz 10 +1287z 8 +924z 6 +330z 4 +45z 2 +1, f 20 (z) = z 19 +18z 17 +136z 15 +560z 13 +1365/ 1 +2002/+ 

1716z 7 + 792 z s +165 z 3 +10 = z(z 2 + 2)(z 4 + 3z 2 + l)(z 4 + 5z 2 + 5)(z 8 + 8z 6 + 19z 4 + 12z 2 +1) . etc . 

Again, we can observe the extensive list above of polynomials in the complex variable. The software 
allows us to list and perform the decomposition (factorization) of high order elements. From an 
appreciation in the list above, we can identify some cases of irreducible polynomials when we have 

the condition F n (z) is ireductible, n = p is prime Qn the other hand, from the list and comparing it 
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with the theorems of the previous section, we can conjecture that the divisibility properties can be 
derived also in the case of the complex variable. Now, from the above relationships, we conjecture the 
following results 

Conjecture 1: For m > 2. we have m,n> 1 F m (z) \ F ;; (z) <^> m\ n . 

Conjetucture 2: For ra, n > 1 then gcd(F„ ( (z), F„ (z)) = F gcd( „, n) (z) . 

Conjetucture 3: F n (z) is irreductible, when n = p is a prime. 


Certainly, the previous conjectures can be verified in some particular cases in order to predict their 
chances of success or the existence of possible counterexamples. This type of procedure or concern is 
recurrent and standard in Mathematics. 



Figura 9. The software allows the determination and decomposition of Fiubonacci polynomial functions in the complex 

variable, (produced by the authors). 

To conclude, in the table below we bring a comparative and conceptual picture of the properties we 
have discussed so far, with emphasis on the Lucas sequence model. In a brief and simplified way, we 
invite the reader to appreciate the existing framework and conceptual relationships closely related to 
the numerical sequences we discussed in the introductory section. We recall the basic properties of 
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1 +V5 rt 1 — y5 2 i r\ 

values a = ———,P =——— the roots of the equation x — x — 1 = 0(Koshy, 2007; 2011; 
Honsberger, 1985; Stakhov, 2005; 2009; Tattersall, 2005). 


Table 1: Comparative table of properties related to the Lucas and Fibonacci model. 


Arithmetic 

properties 

Fibonacci's model 

Lucas' model 

Let k,n e IN , with 
k odd,then 
f \ f, and 

J n J k-n 

L \L, . 

n k-n 

F„ (x, y) \ F (x, y) m\ n 
(Web & Parberry, 1969). 

Corresponding property is not true. 
L 4 (x, y) = x 4 + 4x 2 y + 2y 2 and 
L| 6 (x, y) are irreductible. 

Z^ 9 (x, y) is reductible over Z\x, y ]. 

Let k, n e IN , with 
k even, then 

f n \fkn and 
L n \L k . n . 

F,(x,y)\ F n (x,y) m\n 
(Web & Parberry, 1969). 
F 2 \F <=> F \m 
(Matijasevi'c, 1970). 

Corresponding property is not true. 
A 2 (x, y) = F 4 (x,y)-F 8 (x,y). 

Let m, n e IN and 
let d = gcd (m,ri), 
then 

gcd (f m ,fn) = fd- 

gcd(F„ (x, y), F„ (x, y)) = F gcd(mn) (x, y) 

Corresponding property is not true. 

gcd(L m (x, y), L„ (x, y)) ^ F gcd(m>n) (x, y) 
gcd(L 12 (x, y), L 6 (x, y)) = 1 * L 6 (x, y) 

BLP 

Let m, n e IN and 
let d = gcd(m, n ). If 

m n 

the numbers —, — 
d d 

are both odd, then 
gcd(L m ,L ) = L r 

gcd(F m (x, y), F„ (x, y)) = F gcd(mnj (x, y) 

With the imaginary unit i 2 = — 1. 
BCFP 

Corresponding property is not true. 

gcd(L m (x, y), L„ (x, y)) * L gcd(m>n) (x, y) 

and 

Ig(x,y) = x 8 -8x 6 y + 20x 4 y 2 -16x 2 y 3 +2y 4 

Is irreducitible. 

gcd(L g (x, y), L 4 (x, y)) ^ L 4 (x, y) 

(BCLP) 

_a n -p" 
in ~ a-/3 ’ 
L,=a n +P n 

(a"(x, y)~P"(x, y)) 

F (x, y) = -- - 

a(x,y)~ P(x,y) 

With the imaginary unit i 2 = — 1. 

4 (x, y) = «" (x, y) + P" (x, y) 

With the imaginary unit i 2 = — 1 . 
(BCPL) 

_a n -p n 
in ~ a-p ’ 
L n =a n +P" 

(a" (a, z)-/?"(a, z)) 

F (a, z) = --- 

a(a,z)-/7(a,z) 

Where ‘z’ is a complex variable 
(Taskopru & Altintas, 2015) 

L n (a, z) = a" (a, z) + p" (a, z) 

Where ‘z’ is a complex variable 
(Taskopru & Altintas, 2015) 

f P or L P . where 

p e IN is prime. 

/i 9 =113-37 

F (x, y) and ‘p’ is prime. Then 

F v (x, y) is irreductible over the 
ring Z[x, y]. (Hoggatt & Long, 
1974) 

Corresponding property is not true. 
Z 1 (x,y)\Z 13 (x, y) 

4(x,y)\Z 19 (x,y) 
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(Web & Parberry, 

Conjecture: F (. z) is irreductible 

Lj 2 (x, y) is not divisible by L,(x, y) 

1969). 

over the ring C[x, y] 

and L^x, y). 


3. Conclusion 

In this work we present a set of mathematical definitions related to a generalized sequence model that 
allow an understanding of the evolutionary epistemological and mathematical process of a second 
order recurrent sequence, originally systematically studied by the French mathematician Francois 
Edouard Anatole Lucas (1842 - 1891) and, another formulated by Leornardo Pisano, in 1202. In 

addition, the problem of factorization of numbers present in both sequences {./« }-„ c/v a nd j L n } ^ 

consists of a problem of interest studied in the field of Pure Mathematics and Mathematics Applied for 
decades. Thus, we can find some particular cases in the specialized literature that certain expected 
divisibility properties do not occur, for both numerical sequences (Onphaeng & Pongsriiam, 2014; 
Pongsriiam, 2014). And in many cases, only with the evolution of technology has it become possible 
to test and verify the numerical behavior of a sequence, for increasing indices (see Figure 3). 

Moreover, given these elements and others that we seek to discuss in all sections, mainly some 
elements with respect to an evolutionary epistemological trajectory and, especially, an historical 
perspective (see definitions 1, 2, 3, 4, 5, 6 and 7). In this way, it may raise an understanding about the 
continued ans unstopable progress in Mathematics and some elements, which can contribute to an 
investigation about these sequences, which is customarily discussed in the academic environment, 
however in relation to their formal mathematical value. 

In view of the use of Computational Algebraic System Maple, we have explored particular situations 
enabling a heuristic thought and not completely accurate and precise with respect to certain 
mathematical results (see figures 4, 5, 6, 7 and 8). Such situations involved checking of algebraic 
properties extracted from numerical, algebraic and combinatorial formulations of the Fibonacci’s 
Polynomial model and the Lucas "Polynomials model, for example, discussed by Hoggatt & Bicknell 
(1973a; 1973b). In addition, with the software we can determine the explicit behavior of several 
parficular cases in order to determine possible counterexamples for the both models discussed in this 
work, above all, related to the Lucas' model, for example L,, (x, y) = L 4 (x, y) ■ L%(x, y). 

Another important role is due to the detailed study and epistemological appreciation of the formal 
mathematical definitions addressed throughout the work (set of seven definitions). We assume a 
position that Mathematics progresses and it's progress can be appreciated from the progressive 
establishment of mathematical definitions (Alves, 2017; 2016a; 2016b. 2016c, 2015), since they 
constitute the marked elements of the solidity and certainty of mathematical assertions and therefore of 
theorems discussed here. 

Linally, since all the arguments and properties presented are closely related to the Fibonacci's model, 
an understanding of it's mathematical, epistemological, and evolutionary process cannot be 
disregarded or reduced to a kind of Mathematical History's textbook approach (Arcavi & Isoda, 2007). 
Restricted in an eminently playful and fictional discussion of the Fibonacci episode, in view of the 
problem of the birth of pairs of rabbits. 

In this way, we present to the reader elements and information derived from a formal and 
computational mathematical model, in order to demonstrate a greater, complete and broad 
understanding of some studies on the subject that require a greater scientific discussion and 
dissemination, in order to provide an evolutionary understanding about the current research of some 
subjects in Mathematics. Table 1 should provide a systematic and simplified view of the various 
properties (related to the BFP, BLP, BCFP and BCLP) discussed throughout the paper. 
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